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Abstract 

We reproduce masses of the self-conjugate and non-self-conjugate mesons in 
the context of the spinless Salpeter equation taking into account the rela- 
tivistic kinematics and the quark spins. The hyperfine splittings for the 2S 
charmonium and IS bottomonium are also calculated. Further, the pseu- 
doscalar and vector decay constants of the Be meson and the unperturbed ra- 
dial wave function at the origin are also calculated. We have obtained a local 
equation with a complete relativistic corrections to a class of three attractive 
static interaction potentials of the general form V{r) = —Ar~^ + kt^ + Vq, 
with (3 = 1,1/2, 3/4 decomposed into scalar and vector parts in the form 
Vv{r) = -Ar~^ + (1 - e)Kr'^ and Vs{r) = e/tr^ + Vq; where < e < 1. We 
have used the shifted large-N-expansion technique (SLNET) to solve the re- 
duced equation for the scalar (e = 1), equal mixture of scalar- vector (e = 1/2), 
and vector (e = 0) confinement interaction kernels. The energy eigenvalues 
are carried out up to the third order approximation. 



I. INTRODUCTION 

The spinless Salpeter (SS) equation represents a standard approximation to the Bethe- 
Salpeter equation [1]. Upon elimination of any dependence on time-like variables in a suitable 
manner, the Bethe-Salpeter (BS) equation reduces to the Salpeter equation [2]. Neglecting, 
furthermore, any reference to the spin degrees of freedom and restricting to positive energy 
solutions, one arrives at the SS equation, which is the correct tool to deal with the bound- 
state spectrum of qq, QQ and Qq or qQ (where q = u, d, s and Q = c, b, t) interacting 



via some effective potential with no spin dependence, i.e., spin-averaged data (SAD). This 
equation was solved analytically and then numerically for its bound-state energies using 
different techniques by many authors [3-10]. 

We solved the SS equation [8-10] analytically and then numerically, for the first time, 
by using the shifted large-N-expansion technique (SLNET) for a wide class of static phe- 
nomenological and QCD-motivated potentials previously proposed for quarkonium [11-13]. 
The non-local SS equation was expanded formally to order of f^/c^ up to the first relativistic 
correction term. Thus, the resulting equation is a local Schrodinger-type equation for two 
constituent interacting particles of masses rriq and rriQ. On the other hand, it consists of 
heavy quarks and it can be reliably described by the use of the methods developed for the 
cc and the bb spectra. 

The relativistic effects are not negligible in quarkonia, and we need to investigate them. 
The simplest relativistic effect is just kinematics. They are considered in the potential as 
nonstatic terms, including both spin-dependent and spin-independent parts. Therefore, the 
nonstatic terms in the potential, as in the Fermi-Breit approximation include a spin-spin 
hyperfine interaction, a tensor interaction, a spin-orbit interaction, and a spin-independent 
interaction [14-20]. There are also nonstatic terms in the potential which are independent 
of spin. It is difficult to isolate these terms from data, and furthermore, there is some 
controversy concerning their functional form [21]. In this situation, it is better to omit them 
altogether. 

The discovery of the B^ (the lowest pseudoscalar ^S^ state) was reported in 1998 by the 
collider Detector at Fermilab (CDF) collaboration in 1.8 TeV p-p collisions at the Fermilab 
[22]. The observed mass Mb^ = 6.40 ±0.39 ±0.13 GeV has inspired new theoretical interest 
in the problem [23-31]. Some preliminary estimates of the cb bound state masses have been 
made devoted to the description of the charmonium and bottomonium properties. More 
experimental data are expected to come in near future from new hadron colliders. 

The spectroscopy of the cb system have already been widely studied various times over 
the past years in the framework of the heavy quarkonium theory [23] . The revised analysis 



of the Be spectroscopy has been performed in the framework of the potential approach and 
QCD sum rule [24, 27] . Kwong and Rosner [28] predicted the masses of the lowest cb vector 
(triplet) and pseudoscalar (singlet) states using an empirical mass formula and a logarithmic 
potential. Eichten and Quigg [22] gave a more comprehensive account of the energies and 
properties of the cb system that was based on the QCD-motivated potential of Buchmiiller 
and Tye [32]. Gershtein et al. [29] also published a detailed account of the energies and 
decays of the c6 system using a QCD sum-rule calculation. Baldicchi and Prosperi [27] 
have computed the c6 spectrum based on an effective mass operator with full relativistic 
kinematics. They have also fitted the entire quarkonium spectrum. Fulcher et al. [25] 
extended the treatment of the spin-dependent potentials to the full radiative one-loop level 
and thus included effects of the running coupling constant in these potentials. They have 
also used the renormalization scheme developed by Gupta and Radford [32]. We applied 
the nonrelativistic form of the statistical model [30] and shifted large-N expansion technique 
(SLNET) [31] to calculate the spectroscopy, decay constant and other properties of the heavy 
mesons, including the d) system. 

Recently, in 2002, the ALEPH collaboration has searched for the pseudoscalar bottomo- 
nium meson, the rji, in two-photon interactions at LEP2 with an integrated luminosity of 699 
pb^^ collected at e^e~ centre-of mass energies from 181 GeV to 209 GeV. One candidate 
event is found in the six-charged-particle final state and none in the four-charged-particle 
final state. The candidate rji, {rib -^ KsK^-k^-k^'k^) has reconstructed invariant mass of 
9.30 ±0.02 ±0.02 GeV [33]. Theoretical estimates (from perturbative QCD and lattice non- 
relativistic QCD of the mass splitting between 775(18) and T(1S), M(T(r'^Si)) = 9.460 GeV, 
are reported (cf. [33] and references therein). 

Further, in 2002, the Belle Collaboration [34] has observed a new pseudoscalar char- 
monium state, the 77^(28), in exclusive B — > KKsK~tt~^ decays. The measured mass 
of the r/c(2S), M(r7c(2S)) = 3654 ± 14 MeV. It is close to the r]ci2S) mass observed 
by the same group in the experiment e+e^ — > J/'ipVc where M{ric{2S)) = 3622 ± 12 
MeV was found [35]. It is giving rise to a small hyperfine splitting for the 2S state. 



Ahfs(2S,exp)= M{2^Si) - M{2^So) = 32 ± 14 MeV [36]. Badalian and Bakker [37] cal- 
culated the hyperfine splitting for the 2S charmonium state in a recent work. Recksiegel 
and Sumino developed a new formalism [38] based on perturbative QCD to compute the 
hyperfine splittings of the bottomonium spectrum as well as the fine and hyperfine splittings 
of the charmonium spectrum [39] . 

The motivation of the present calculation is to extend the SLNET [8-10,31,40-42] to the 
treatment of the SS equation [10] to calculate the hyperfine splittings for the 2S charmonium 
and IS bottomonium and also to reproduce the cb Salpeter binding masses below the con- 
tinuum threshold M{AS). We shall use a potential model that includes a running coupling 
constant effects in the central potential to give a simultaneous account of the properties of 
the cb system. We choose a class of three static central potentials [10-13,30,31,41,42] each 
with a strong coupling constant a^ to fit the spectroscopy of the existing and non-existing 
quarkonium systems. Since one would expect the average values of the momentum trans- 
fer in the various quark-antiquark states to be different, some variation in the values of 
the strong coupling constant and the normalization scale in the spin-dependent should be 
expected. We extend our previous work [31] for the solution of the SS equation to deter- 
mine the binding masses of the cc, bb, and cb mesons by taking into account the spin-spin, 
spin-orbit and tensor interactions [14-21]. The spin effects are treated, in the framework of 
perturbation theory, as perturbation to the static potential and the treatments are based on 
the reduced Salpeter equation. 

The outline of this paper is as following: In Section II, we first review briefiy the analytic 
solution of the SS equation, we consider here the unequal mass case studied in our earlier 
works for a class of static potentials using SLNET. Section III is devoted for the class of three 
static potentials, which are decomposed into scalar and vector parts and also for their spin 
corrections. Pure vector and scalar potentials as well as their mixture are considered. Section 
IV is devoted for the pseudoscalar and vector leptonic constant of the Be meson. We finally 
present an approximation to cc, bb, and cb mass spectra and decay constant of the Be meson. 
We also calculate the the observed new charmonium ric{2S) and searched bottomonium 



?7c(lS) mesons and the hyperfine splittings of their states. Finally, our conclusions are given 
in Section V. Appendix A contains some definitions as well as the formulas necessary to 
carry out the above mentioned computations. 

II. WAVE EQUATION 

The relativistic wave Salpeter equation [8-10] is constructed by considering the kinetic 
energies of the constituents and the interaction potential. We assume that the full Hamilto- 
nian H governing the dynamics of the quantum-mechanical system under consideration can 
be split up into a free Hamiltonian Hq and an interaction potential V{r). For the case of two 
particles with unequal masses iriq and mg, interacting via a spherically symmetric potential 
y(r), the spatial coordinate representation of the SS equation in the center-of-momentum 
system for the two-body system reads off 



'<-=q,Q 



m (r) = 0, (1) 



with the potential V{r) contains, in addition to the static interaction Vstatic{f^) the total 
spin-dependent potential Vsoif^) [15-18,20,21]. The kinetic terms involving the operation 



Ajv + iTif are nonlocal operators and can be solved through a direct expansion of the 
square root-operator [5,10,17]. For heavy quarks, the kinetic energy operators in Eq. (1) 
can be approximated by the standard expansion in the inverse powers of quark masses to 
obtain the usual nonrelativistic Schrodinger Hamiltonian [5,10] 



Aa^ a 



2 



E ^-^N + m^ = m, + mQ-^-^^^...^ (2) 



i=q 



where u = "^'^"^'^ denotes the reduced mass and ri = a ( — ^i"^Q ] jg a useful mass 



,1/3 
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parameter. The radial part of Eq. (1), in the N-dimensional space, is expanded in powers 
of f^/c^ up to one relativistic correction term as [5,8-10,16] 



Am a 



N ^N 



3 



2/x 8r] 



\/(r) i?„,i(r) = K,;i?„XO, (3) 



where En,i = Mn,i{qQ) — rriq — itlq stands for the Salpeter binding energy and A 



N 



N- 



^ This SS-type equation retains its relativistic kinematics and is suitable for describing the 
spin-averaged spectrum of two bound quarks of masses rriq and rriQ and total binding meson 
mass eigenstate Mn^i{qQ). Furthermore, in order to obtain a Schrodinger-like equation, the 
perturbed term in Eq. (3) is treated using the reduced Schrodinger equation [44] 



p^ = 4/i2 [E^^z - V{r)f , 



(4) 



where p is the center of mass momentum of the qQ system. Consequently one would reduce 
Eq. (3) to the Schrodinger- type form [17] 



[ 2/i 2r]^ 
After employing the following transformation 



Ki + ^ (^) - '^En,iVir) + Vir) i?„,Kr) = ^„,/^n/r) 



(5) 



RnAr) 



r(^-i)/2' 



(6) 



we may rewrite Eq.(5) simply as 

2fidr'^ 
with 



'k-{l-a) 




'k-{3- a) 



i^r^ 



+ Wr^jir) - 



Wn,i{r) 
2m 



Un,i{r) = 0, 



(7) 



Wn,i{r) = Veffir)-En,i, 



(8) 



and the effective mass parameter 



m =r] /fi = {mqmQjj,)/{mqmQ - 3/i ). 



(9) 



The perturbation term, Wn^iir^, in (7) is significant only where it is small (i.e., Wn^i{r)/m ^ 
1). This condition is verified by the confining potentials used to describe heavy— quark 



^This approximation is correct to 0{v'^/c'^). The A^ term in (3) should be properly treated as a 
perturbation by using trial wavefunctions [43] . 
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systems except near the color— Coulomb singularity at the origin, and for r -^ oo. However, 
it is always satisfied on the average as stated by Durand et al. [17]. We now proceed to 
solve Eq. (7) with Eq. (8). In the SLNET [8-10,30,40-42], it is convenient to shift the origin 
of coordinates to r = ro (or y = 0) hj defining 



y = k^/'^{r - ro)/ro. 



(10) 



Expansions about this point in powers of y and k yield 



1 - i f ,_,).(i±i) 



r^ r. 



j=o 



k 



-i/2 



(11) 



V{y{ro)) 



1 ^ ( d^V{ro) \ (ro^y (4-,)/2 



(12) 



and also 



E,^ 



nl 



1 ,^^2-i 



Further, by substituting Eqs. (8) through (13) into Eq. (7), one gets 



(13) 



1 (P 



4:fidy' 



+ 



k (2 -a) (l-a)(3-a) 



16/i 8/x 



16fik 



^f ,,(1+J) 

X L i-y) ■ 



j=0 



rj/2 



Q h V dri ) j\ m'Q 



- ( d^V{r,) \ MV p-,/2 



2rl - (1-,) - f dW{ro) \ {royU-j/2 . 



(14) 



with the eigenvalues 



E' 



Q [ \ mix m 



2EqEi 



„ 2EnEo E-i \ 1 

E^ + ^^ + ^ r 

m m k 



, „ lEnE'^ 2E-\ E<2 \ 1 
m m J k'^ 



(15) 



The parameter Q is an arbitrary scale, but is to be set equal k at the end of our calculations. 
Thus comparing Eqs. (14) and (15) with its counterpart Schrodinger— like equation for the 
one— dimensional anharmonic oscillator problem [40], we calculate all the relevant quantities 
(5's and e's [31]. The final analytic expression in the Xjk expansion of the energy eigenvalues 
appropriate to the SS particle is 



^Ur ~ k 



1 rgl^(ro) rlV{r^f 2rlE^V{r^) 



16/i 



Q 



m'Q 



m'Q 



+ 



:i + 2nr 



UJ 



(2-a) 
8/x 



IrlE^Vjr^) 
mQ 



+ /9 



(1) 



1 



2rlE,V{r,) 
m'Q 



+ /9 



(2) 



+ 



1 
F 



(16) 



with rij. is to be set equal to 0, 1, 2, ■ ■ ■ . The quantities 13^^^ and /S^^^ appearing in the correc- 
tion to the leading order of the energy expression are displayed in [31]. 

Comparing the terms of Eq. (15) with their counterparts in Eq. (16) and equating terms 
of same order in /c, one gets the leading contribution to the binding energy 



m 



Eo = V{ro) + — 



'1 + 



Q 



ArQin^m"^ 



-1 



(17) 



where m^=/i/'m . Here r^ is chosen to minimize Eq, that is. 



dEo 
dro 



d^Eo 

"■• id"'' 



(18) 



and it satisfies 



roV'(ro) 



m 



Q 



1/2 



Q 



IQrnin 



0'"-^. 



16m, 



(19) 



To solve for the shifting parameter a, the next contribution to the energy eigenvalue, in Eq. 
(13), is chosen to vanish (i.e., Ei = 0) which gives 



a = 2- 4/i(l + 2nr)uj, 



(20) 



with UJ defined as 



UJ 



4/i 



3 + roV^"(ro)/V^'(ro) - \^rlm^V'{r^f jQ 



1/2 



(21) 



The scaling parameter Q, in Eq. (19), is simply written as 



Q = 8m^ roV(ro) (1 + 0, 



with 



(22) 



i 



m 



(23) 



Therefore, with the help of relations (20) through (23) together with Q = k^ , we obtain the 
following formula 



1-2 



1 + 21 + 4/i(2n, + l)uj = 2r^\/'(ro) (2m^ + 2m^0 



,1/2 



(24) 



which is an explicit equation in tq. Once r^ is determined via (24), it becomes easy and 
straightforward to obtain Eq via Eq. (17), E2 and E^ via solving Eqs. (15) -(16). Thus, 
the general expression for the binding energy eigenvalue takes 

1 



Er, J = En 



rin- 



•2W„,liro) 



'iP'^SL^ol' 



k 



A;2 



(25) 



which works well and is convergent as the value of / increases. Finally, the binding meson 
mass eigenstate for the quarkonium families is 



Mn,l{qQ) = '^En,l + mq + mQ. 



where mi and 7712 are constituent quark masses. 



III. HEAVY QUARKONIUM AND Be MESON MASS SPECTRA 



(26) 



To describe the spin-dependent relativistic corrections to the potential, we have used a 
similar approach derived by Olsson et al. [14] for positron. We choose the potential in Eq. 

(1) as [18,20,44] 



V{r) = Vstaucir) + VsD{r) + Vsi{r), 



(27) 
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with spin-dependent and spin-independent perturbation terms are given in Refs. 
[20,21,43,44]. Further, the static potential [7-13,30,31,41,42] takes the general form 

Vstatic{r) = -—p+ >^r^ + Vo; 13 = 1, 1/2, 3/4, (28) 

where A > 0, /t > and Vq "^^J be of either sign. The first potential we consider here is 
the Cornell [13] potential which is one of the earliest QCD-motivated potentials having the 
form 

Vc{r) = -^ + Kr + V^, (29) 

where A = 4as/3, the Coulomb-type piece characterizes the short-range gluon exchange, 
K is a confinement constant, and the constant Vq is for spin-independent interactions not 
included explicitly in the r— dependent part of the potential. It is related to the slope k, of 
the linear potential by [15] 

K) = -2v^exp(-7s + -), (30) 

where 'Je = 0.577215 ■ ■ • is the Euler-Mascheroni constant. The second potential is that of 
Song and Lin [12] and is given by 

Vs-L{r) = -^ + Kri/2 ^ y^_ (31) 

The third potential is an intermediate case between the last two mentioned potentials and 
is called Turin potential [11] 

VT{r) = -^^+Kr'/' + Vo. (32) 

The class of static potentials in Eq. (28) must satisfy the following conditions 



dV cPV 

dr ' dr 



— > 0' ^ < 0- (33) 



The fine and hyperfine splittings are computed by using the Breit-Fermi interaction with 
the assumptionp [16,18] 

11 



Vs(r) = enr^ + Vc 



(34) 



Vvir) = VstaUcir)-Vsir), 



(35) 



where e is some mixing parameter. The vector term incorporates the expected short-distance 
behaviour from single-gluon exchange. We have also included a multiple of the long-range 
interaction in Vv{r) to see if we can determine the vector-scalar nature of the confining 
interaction. Here we investigate the cases of pure scalar confinement (e = 1), equal mixture 
of scalar- vector couplings (e = 1/2) and a pure vector case (e = 0). The spin dependent 
correction to the nonrelativistic Hamiltonian, which is responsible for the fine splittings, is 
also modelled on the Breit-Fermi Hamiltonian [22,29]. It can be decomposed into a part, 
which is antisymmetric with respect to the spins of the constituents Vyi(r), and apart sym- 
metric in these spins. The symmetric part can be decomposed into a spin-orbit interaction 
VLsif^), ^ spin-spin interaction Vss{i^), and a tensor part Vrir). Therefore, it is given by 
[15-17,20,21,44] 



VsDir) = VA + Vs 



mi 



7717 



Vl.{r) - V^{r) 



L S 



1 L-S V{,{r) 1 1 


L ■ Si L ■ S2 

_ m2 ml _ 




V^{r)-V',{r) 


mqm,Q r 2 


r 


2Si-S2 
Smgmg 


VVy(r)' 


S12 
nigmq 


-Vy 


(^ 


-, . ^^c-)] 


1 



(36) 

where Si and S2 are the quark spins, S_ = Si — S2, L = x x p is the relative orbital angular 
momentum, and S12 = T — (Si ■ S2) /3 where T = (Si ■ f)(S2 ■ f) is the tensor operator 
with the versor r = r/r . The spin dependent correction (36), which is responsible for the 
hyperfine splitting of the mass levels, in the short-range is generally used in the form for 

^-wave (L = 0) (cf. e.g., [32,43,44]) 



VHFs{r) = -{Si-S2)V' 



4a. 



3^/3 



(37) 
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but the one responsible for the fine sphttings is used for P- and D-waves {L ^ 0) it is given 
by 

1 fL-S 



VFs{r) 



rriqmQ { r 



1 + -.— I W/fr) - -^ ^V^^frl 



4 m^mg 



V 



4 mgrng 



+ 3(81 ■S2)V^ 



kCI — t)r- 



T--{S,-S,] 



K'(r) + 



\/^(r) 



(38) 



where the matrix element can be evaluated in terms of the expectation values [22] (L ■ Si) = 
(L ■ S2) = I (L ■ S). Hence, Eq. (38) is the complete spin-dependent potential in QCD 
through order m^. For bound state constituents of spin Si = S2 = 1/2, the scalar product 
of their spins Si ■ S2 and L ■ S are to be found in the Appendix A. The appearance of a 
Coulomb-like contribution ~ 1/r in the vector part of the potential causes some problems 
due to the relation V^(l/r) = —A7r5^^\x), in the spin-spin interaction (37) involves a delta 
function of the S'-wave (L = 0). Thus, for Cornell potential, the hyperfine splitting potential 
(37) gives 



V,ff{r) = + Kr + -^^^5(3) (r) (Si ■ S2) + Vo; where (3 = 1. 

r drnqiriQ 

Therefore for the energy of spin-spin interaction we have approximately: 

1 



(39) 



E., 



2M, 



n,0 



-AMi (Si • S2) , 



(40) 



where M„ is given in Eq. (26) and the singlet-triplet mass squared difference [15] 



AM^ 



32 



(41) 



for light qq systems (in the instantaneous-limit approximation) [15], and 



am! 



m: 



S=l 



, .2 256 



(42) 



for heavy quarkonia (hydrogen- like trial functions) [15]. All these predictions for the mass- 
squared difference are independent of the mass of the particles which constitute the bound 
state. Further, for the Song-Lin and Turin potentials, it also give 
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V,ff{T) = -4 + >^r^ + ML^lZ^r-'^Si ■ S2 + 1^0; where (3 = 1/2, 3/4. (43) 

On the other hand, the Eq. (38), for P, D,- ■ ■ waves (L 7^ 0) case, gives 

Veffir) = Vstaucir) + 9{r) \Fls_ (L ■ S_) + Fls (L ■ S) + F55 (Si ■ S2) + FtTI , (44) 



with a given set of spin-dependent quantities 



'LS- 



4 rrigmQ 



Ar 



-13 



'1 — e) nr 



(45) 



'LS 



1 + 



1ml + ml' 
4 m^mQ ^ 



Ar~ 



:i - 2e) Kr 



e/tr'^ 



(46) 



and 



'ss 



(2 + /9) 



Ar"'^ + /3 (1 - e) Kr/^ 



(2 + /5)Ar~^ + (2 - /3) (1 - e) ^r^' 



(47) 



(48) 



where g{r) 



mamnr 



.qmqr 



is a necessary couphng function. The spin-independent corrections in 



Eq. (27) are exphcitly given by Ref. [44] which are not treated in our present work. 



A. Singlet states 

For parastates {L = J) or {S = 0) case, we have parity P = (— 1)^^+^ and charge conju- 
gation C = (—1)"^. Thus, the potential (44) can be rewritten as 



Veffir) = Vstaticir) - - {3Fss + Ft) 



(2L + 3)(2L-l)Fi5_ 



10 



(49) 



which can be substituted in Eq. (8) and also by setting k = N + 2J — a therein. Further, 
Eqs. (39) and (43) give 



Veffir) = -- + Kr 



87rQ;s 
Smgmg 



S^'\r) + Vo, 



(50) 
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and 



Veffir) = -4 + -^ - '^ii/jr^ ^'' + ^0' -here /5 = 1/2, 3/4, 



r./3 



(51) 



respectively, which generate singlet states with opposite quark and antiquark spins of the 
signature u^Sq. Furthermore, Eq. (49) can be rewritten simply as 



Vj=L{r) = g(r) < - ^ 



4 mqmq V 10 



— (2L + 3)(2L - 1) \Ar-l^ + (1 - e)/€r^ 



-Xl + P){l-e)KA + VstaUc{r) 



(52) 



which generates states of the signatures n^Pi, n^D2, n^F^, ri^Gi,---. 



B. Triplet states 

For triplet {S = 1) case, we have the known inequality \L — S\<J<L + S that gives 
J = L and L ± 1 : 

i. States J = L 



Here, the parity P = (— l)"'+^and the charge conjugation C = (—1)^^^. The potential in 
Eq. (44) takes the following simple form 



Veffir) = Vstaucir) + \ {Fss + Ft - AFls) + \l ^{21 + 3)(2L - 1)Fls_: 



(53) 



4 ' ^^. y ^Q 

which can be substituted in (8) together with k = N + 2J — a therein. Further, the potential 
(53) reads 

g{r 



Vj=Lir) 



Iml + ml 



4^ l^ll^U.- /, 
^ 3 2 mgiriQ J \ 2 rUgmQ 



'1 — 2e)Kr^ + enr 



.P 



1 m% — rn? 



-dj^C^L + 3)(2L - 1) [Ar-f^ + (1 - e)Kr 



2 mqmQ 
which generates states like n^Pi, n^D2, n^F^, n^G^ 



+ Vs 



static) 



(54) 
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2. States J = L±l 

We have the parity P = (—1)"^ and the charge conjugation C = (—1)^+^. The eigenfunc- 
tion is a superposition of two components with orbital momentum L = J+ 1 and L = J — 1 
which have equal space parity 



^sAr) = Uj.^{r)Yy^_,^,A^,v)+Uj+^{r)Yj"l,^,A(^,v). 



(55) 



The action of the tensor operator, T, on the two components of the wavefunction in Eq. 
(55) is 



Tuj^iyr±i,i,ji^) = Ti(^7^«^±i^i^i,w(?) + l \j+i^ '''^^'^'^Ti,U'^)- (56) 



Therefore, a set of equations are obtained 
2/i dr'^ ' Syur^ 



'k-{l-a) 




'k-{3- a) 



1 / Ft 

Vstatic{r) + - I Fss - /2j _L. i) 



E. 



n,J+l 



(J + 2)F..-- 



Vstaucir) - ( J + 2) Fls + ^ [Fss - (^J^ 



E^, 



n,J+l 



IJ{J+l) 



(57) 



and 



1 (f 



2/i dr' 



+ 



'k-{l-a) 




'k-{3- a) 



8/ir^ 



1 / -Ft \ 



+ (^_i,f,,__ 



V;,<.,.o(r) + (J - 1) Fis + i (fss + pj^ 



H 2 



En, J -I 



IJ{J+l) 
XUn,J-l[r) - ' FTUn,.J+l{r) = 0, 



(58) 



2(2J+1) 

where A; = A^ + 2J + 2 — a. Therefore, Eqs. (57) and (58) describe states such as 
n^P2, n^D^, n^F2, n^H4^,n^PQ, n^Di,---. Here we may consider numerically the system 
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obtained and separate equations by dropping out the mixed terms to see their effect on 
the spectrum of the masses. Consequently one can rewrite (57) and (58) in the following 
simplest forms 



Vj=L-iir) = -gir)\iL + l) 



\^}_^l + ^Q 



4 mgrriQ 



Ar^^ + (1 - 2e) kt^' 



eKV 



+- 



1 1 



4 (2L - 1^ 



(2 + f3)Ar-'^ + (2 - /3) (1 - e) /tr^' 



1 

+ 4 



(2 + /5) 



Ar-^ - l3 {1 - e) Krf^ 



+ Vstaticif 



staticyi ) i 



for states n^Po? n^Di, n^F2, n^H^, ■ ■ ■ and 
Vj=L+i{r) = g{r)- 



l + i!"9 + "^Q^ 



4 rnqruQ 



Ar'^ + (1 - 2e) nr^ 



+ enr 



L 



(59) 



+ - 



1 1 



4 (2L + 3) 



(2 + /?)Ar-'^ + (2 - /?) (1 - e) nr^ 



(2 + /5) 



Ar-'^-/5(l-e)Kr'^ 



+ Vstatic{r), 



for states n^P2, n^D^, ■ ■ ■ . Further, for triplet S'-wave, we have 

Veffir) = + Kr + ^^^^'^ r + Vo, 

r 9mgmQ 



and 



Veffir) 



A 



p 2/J(l-/J)7ra,^.^ 



kt'' + 



which describe states such as n^Si, 



r'" + Vo; where /5 = 1/2, 3/4, 



(60) 



(61) 



(62) 



5. State J = 



Equations (59) and (60) degenerate into a single equation with an effective potential 



Veffir) = Vstaucir) + -{Fss -Ft- 8Fls), 



(63) 
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and also by setting k = N + 2 — a therein. Further, Eq. (63) becomes 

,2 I ^2 



+ 1 



Vj=o(r) 



gyr) 



1 +g{r) 



2^1S + "^Q , 2 + /3 



2 mqUiQ 



+ 



Ar" 



,5-^)(l-.)+f<t^Vl-2.) 



TUgTUQ 



KT 



' + v,. 



(64) 



which only describes states such as tt'Pq. 



IV. PSEUDOSCALAR AND VECTOR DECAY CONSTANTS OF THE Bq 

MESON 

The significant contribution to the B^. total decay rate comes from the annihilation of 
the c quark and h antiquark into the vector boson W^ which decays into a lepton and a 
neutrino or a quark-antiquark pair. 

The nonrelativistic expression for the decay constants is given by [45-47] 



/. 



NR _ rNR 

— Jv 
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1^ 



py 



(65) 



V MpyiqQ) 

where \E'p,y(0) is the meson wave function at the origin r = 0. The fp and fy, P corresponds 
to the pseudoscalar Be and V to to the vector _B* mesons and Mpy{qQ) are the masses of 



the Be and -B* mesons. 



V. RESULTS AND CONCLUSION 



We have solved the spinless Salpeter equation using SLNET and also extended our ear- 
lier formalism for the SAD spectra [31] by introducing the spin corrections. Furthermore, 
we have obtained a unified description of the self-conjugate meson spectroscopy with a phe- 
nomenological and a QCD-motivated potential model. This model has also extended to 
comprise various cases of pure scalar confinement (e = 1), scalar- vector couplings (e = 1/2) 
and the vector confinement (e = 0) interactions. It was motivated by our desire to construct 



an analytic expression for the quarkonium mass spectra of different spins and confining 
interactions that could have combined both heavy and light quarkonia. 

The parameter fits used in this work together with the quark masses are shown in Table 
I. They are taken to be same as those considered in our earlier works [10] for each static 
attractive potential. We do not view our present result as a totally successful model of 
quarkonia interactions even though the level structure is good and the transitions are also 
well accounted for. It seems likely that in this study our model has shown remarkable success 
in describing the heavy quarkonia data well once the potential parameters are fitted properly 
for optimum agreement with the experimentally observed data. 

We have made the Salpeter binding mass Eq. (26) more general by including the spin- 
dependent corrections to the potential. Such a formula is able to describe with good accuracy 
the spectra of all quark- ant iquark bound states as long as the mass of the constituent 
quark is larger than the term Wn,i{r) as remarked earlier. Our predictions for the Salpeter 
mass spectrum of all unobserved and observed charmonium and upsilon systems, to all 
static potentials, in the flavour-dependent and flavour-independent cases with e = 0, 1/2, 1 
are presented in Tables II-III. Further, Table IV shows the predicted Salpeter masses of 
unobserved n^^^^Dj, ri^^^^Fj and n^^^^Gj bottomonium levels using a new set of potential 
parameters. Therefore the calculated bottom masses of all static potentials for the n^^^^Sj, 
fi^s+ipj and n^^^^Dj states are quite fair for the Turin potential as clearly seen from 
the average relative error range, Um"^^ — M^^'^) / M^'^^'^y We have also predicted the Be 
meson spectrum for the three static potentials in Table V using the set of parameters in 
Table I. In Tables VI and VII we calculated the c6, cc and hb using a different set of fitted 
parameters for Cornell potential. No additive constant is permitted in this manner and the 
accuracy in producing spectra is pretty good. There is a clear preference in these fits for the 
Salpeter wave equation with an approximately equal mixture of scalar and vector (e = 1/2) 
couplings and also vector confinement (e = 0). We found that larger values of e are mildly 
preferred by fits to the spin-singlet or triplet states while the quality of the e = 1 fit to 
the spin singlet and triplet data is not as high as that of the e = 1/2 and e = fits, it is 
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certainly acceptable, with a value of errors. We conclude here that the Lorentz structure 
of the confining interaction cannot be determined using only e = 1/2. We will henceforth 
restrict our attention to the case of pure scalar confinement at large distance as expected 
theoretically. The three potentials seem to be fairly good in fitting all the data.. 

In the equal scalar and vector couplings, we have found that our fits are very good with 
level values and accurate to a few MeV. For convenience we compare explicitly the predicted 
and measured spin splitting energy for different L states. We find that the apparent success 
is achieved for the predicted Xb2 — Xbi = 24 MeV and Xbi — Xbo = 33 MeV in the average 
for the three potentials and are very close to the experimental values 21 MeV and 32 MeV 
respectively. Furthermore, the predicted x'b2 '~ x'bi — 13 MeV and x'ti ~ x'bo ~ ^^ MeV in the 
average for the three potentials which are exactly same as the experimental value 13 MeV 
and close to 23 MeV, respectively. The predicted hyperfine splitting Ahfs(lS) = M(T(1S)) — 
M(r/b(lS)) = 80tt MeV, (cf. [33]), Ahfs(2S) = M(T'(2S)) - M(r?^2S)) = 2211 MeV, and 
Ahfs(3S) = M(T"(3S)) - M(r/^' (3S)) = Utl MeV are nearly close to the theoretically 
calculated values 62 MeV, 40 MeV, and 15 MeV, respectively, (cf. [33,39]). Further, The 
hyperfine splitting for the 2S charmonium state is calculated and the predicted number is 
Ahfs(2S) = M(t/'(2S))-M(r?c(2S)) = 56l^ MeV for fiavour dependent case and 56t^^ MeV for 
fiavour independent case, (cf. [37,39]) . Badalian and Bakker in their recent work calculated 
and predicted the number as Ahfs(2S, theory) =57 ± 8 MeV [37] giving M{ric{2S) =3630±8 
MeV. 

We observe, in this connection, that splitting approximation can be improved signifi- 
cantly by increasing the quantum number L. It has also been noted that the term contribut- 
ing most to the Salpeter binding energy is the leading term Eq as L becomes non zero which 
improves the convergence of Eq. (25). This is expected since the expansion parameter 1/k 
becomes smaller as L becomes larger. 

It is found that most results agree moderately well with the experimental data, or with 
the theoretical predictions wherever experimental data are not available. Our predictions 
for the lower states are considerably better than the higher states since T(nS), where n > 4 
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lie above the threshold and predominantly decay into mesons containing charm and beauty 
flavour.^ 

The results obtained are a relevant summary of the effectiveness of the method. It 
seems likely, that the average relative errors of the fitting results are almost very small for 
most states. The comparisons with the experimental data show clearly that when we are 
taking into account the relativistic corrections and also the third order approximation to 
the binding energy, the coincidence with the experimental data is improved. The deviations 
from experiment are more considerable at n^So states for cc and bb systems. These results 
show that the last terms in these equations in most cases are not important for the spectrum 
of meson masses. We do not pretend to display a whole account of the meson spectroscopy. 
Our purpose is to demonstrate the possibility of applicance of SLNET for determination of 
meson masses in the context of the SS equation using relativistic kinematics. 

We have found that fitting the parameters is extremely essential to enable one to sharpen 
the analysis.'^ The prescription and the wave equation used are responsible for the deterio- 
ration of the fit parameters. The calculation and parameters are also model dependent as 
remarked in our earlier papers [10,31]. All we can say that the SLNET works for the SS 
equation well with relativistic kinematics as these parameters are fitted properly and once 
the mass of the quark is taken large. 

^From a phenomenological point of view, a certain amount of flavour dependence is 
required. It was pointed out by Miller et al. in [14] that no flavour-independent potential 
could fit both the spin-averaged cc and hb levels. 

It is clear that the coulomb-like parameter A is in accordance with the ideas of asymptotic 
freedom is expected for the strong gauge-couphng constant of QCD, that is, asirric) 7^ 



^In the upsilon family, Kiselev [48] found n = 4. 

^We have used the available parameter fits in the literature as a test for the present work. 
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ttsl^fc) 7^ ttsl/u) with /i = merrily/ {rric + mt); (cf. e.g., Ref. [48]).^ Of course a complete 
analysis would require ingredients that we have not considered here but discussed in the 
literature. 

The calculated values of the pseudoscalar and vector decay constants of the B^ meson 
in our model using the nonrelativistic formula (65) are displayed in Table VIll. The radial 
wave function at the origin is being calculated in Table IX. They are compared with the 
ones calculated using the relativistic and other predictions using the nonrelativistic quark 
models. The calculated values of these leptonic constants are consistent with the other pre- 
dictions [22,25,29,38,49,50]. Kiselev et al. [49] estimated the pseudoscalar leptonic constant 
f^^ = 493 MeV in the potential model. Furthermore, the calculations in the same potential 
model with the one-loop matching [50] is f^~ °°^ = 400 ± 45 MeV and also for the two-loop 
calculations [50] is fl'J""^ = 395 ± 15 MeV. 

We may expect to improve our results by adding extra spin-independent terms compa- 
rable to the ones induced by the original scalar confinement potential. 
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^For a pretty good fit to the cc, bb and cb quarkonium spectra, the QCD coupling constant asifJp) 
must be dependent on the quark- flavour. Kiselev [48] took as{m{,)/as{mc) — 3/4. 
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APPENDIX A: THE SPIN-CORRECTION TERMS: 



For parastates {S = 0) case we have: 



J = L 



(Al) 



For triplet (S* = 1) case we have the following: 

L-1, S-L = -(L + 1) 
J=i L, S-L = -l (A2) 

L + 1, S-L = L 

For bound-state constituents of spin Si = S2 = 1/2, the independent operators Si • 
S2,:(Si ± S2) ■ L and T : 



\Si ■ S2) 



-3/4 for spin singlets, S* = 0, 
-1/4 for spin triplets, S* = 1. 



(A3) 



(S-L) 



0, for spin singlets S = 0, 

\ [J{J H-l) - L{L + 1) - 2] , for spin triplets S = 1. 



(A4) 



[Si ■fS2-r)Mj(r)Y'7o,j(f) 



--Mj(r)Fj,o,j(f) 



(Si ■ S2)lX^,Jf) = - [S{S + 1) - Si{Si -f 1) - ^2(^2 + 1)] lX^,Jr), 



[Si + S2) ■ L1X5,l(?) = 2 y^J + 1) - L{L + 1) - S{S + 1)] y}- , Jr) 



(A5) 
(A6) 
(AT) 



Si - S2) ■ ^Y1:s,l{^)=^Yq [2^ + 3)(21^ - l)]5j,L (5s,olXi,L(f) + 5s,iY^,A^)) , (A8) 



(L + l)(L + 2) \L{L-l) 

f" '\/'Tfi /-^A V f" ~\/"tTi /-^A 

2(2L + 3) '^■^'^+1^-^.1-^+2^ - 2(2L-1)'^-''^-' ■^'i'^-^^'"^- 



(A9) 
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TABLES 



TABLE L Fitted parameters of the class of static central potentials. 



Cornell 



Song-Lin 



Turin 



cc 

rfic 

A 

K 



L84 GeV 
0.509 

0.166 GeV^ 
-0.799 GeV 



L80 GeV 
0.498 Gey V2 
0.694 Gey3/2 
-L314 GeV 



L82 GeV 
0.547 GeV^/^ 
0.309 GeV'^/^ 
-0.898 GeV 



bb 

rrib 

A 

K 

Vo 



5.17 GeV 
0.443 

0.184 GeV^ 
-0.805 GeV 



5.20 GeV 
0.762 Geyi/2 
0.576 Gey3/2 
-0.991 Gey 



5.18 GeF 
0.576 GeF 1/4 

0.304 Gey^/^ 
-0.845 GeV 



Combined 

rric 

rub 

A 

K 



1.79 GeV 
5.17 GeV 
0.457 

0.182 GeF^ 
-0.790 GeV 



1.82 Gey 
5.20 Gey 
0.722 Geyi/2 
0.589 Gey3/2 
-1.029 GeV 



1.80 Gey 
5.18 Gey 
0.578 GeV^/^ 
0.304 Gey^/^ 
-0.838 GeV 
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TABLE II. Salpeter masses of the unobserved cc and hh for the flavor-dependent and fla- 
vor-independent cases, respectively for the Cornell, Song-Lin and Turin potentials 



State^ 


[45] 


Cornell'^ 


Song-Lin 


Turin 


Cornell 


Song-Lin'^ 


Turin'^ 


r/c(25)d 


3589 


3604 


3595 


3593 


3602 


3574 


3585 


{l^D2)cc 


3826 


3755 


3756 


3756 


3748 


3756 


3754 


{l^Di)cc 


3814 


3708 


3708 


3705 


3699 


3703 


3701 


%(15) 


9314 


9398 


9407 


9395 


9380 


9433 


9399 


%(25) 


9931 


10000 


10000 


9997 


10000 


10014 


10003 


r]b{2,S) 


10288 


10342 


10329 


10332 


10343 


10341 


10338 


Vhi^S) 


10577 


10612 


10570 


10586 


10613 


10581 


10592 


(liPi)66 


9906 


9905 


9890 


9897 


9908 


9901 


9902 


(2iPi)66 


10254 


10246 


10246 


10245 


10249 


10256 


10251 


{l^D2)hh 


10155 


10142 


10149 


10146 


10145 


10157 


10152 


{2^D2)hh 


10450 


10424 


10422 


10423 


10427 


10431 


10429 


X6o(3P) 


10519 


10512 


10491 


10498 


10513 


10501 


10504 


X6i(3P) 


10542 


10520 


10500 


10508 


10522 


10510 


10514 


X62(3P) 


10561 


10528 


10508 


10516 


10530 


10518 


10522 


{l^Di)bh 


10147 


10128 


10132 


10131 


10131 


10142 


10137 


{l^D2)bb 


10153 


10140 


10147 


10145 


10144 


10156 


10151 


{l'D3)bb 


10158 


10152 


10161 


10158 


10156 


10170 


10164 


{2^Di)bb 


10442 


10414 


10411 


10412 


10417 


10420 


10418 


(2^D2)bb 


10448 


10423 


10421 


10422 


10426 


10430 


10428 


{2^D3)bb 


10453 


10433 


10431 


10432 


10436 


10440 


10438 



'^States with equal mixtures of scalar and vector couplings. 

Here's the flavour-dependent case. 
*^Here's the flavour-independent case. 
'^Recently observed by Belle Collaboration [34,35]. 
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TABLE III. Salpeter cc and 66 quarkonium masses of the observed states (in MeV) for the 
Cornell, Song-Lin and Turin potentials 



State^ 


M(^^p)[36] 


[45] 


Cornell'^ 


Song-Lin 


Turin^ 


CornelF 


Song-Lin^ 


Turin'^ 


VcilS) 


2980 


2967 


2981 


2965 


2950 


2988 


2983^^ 


2999^^ 


V'(i5) 


3097 


3126 


3104 


3108 


3111 


3092 


3096^^ 


3087^^ 


V'(25) 


3686 


3700 


3652 


3655 


3652 


3650 


3640 


3639 


XcoilP) 


3415 


3427 


3316 


3374 


3335 


3322 


3340 


3309 


XciilP) 


3511 


3497 


3471 


3461 


3464 


3456 


3463 


3459 


Xc2{lP) 


3556 


3543 


3535 


3521 


3530 


3521 


3532 


3529 


hcilP) 


3526 


3510 


3488 


3468 


3478 


3471 


3475 


3474 


TilS) 


9460 


9511 


9474 


9499 


9488 


9466 


9505^^ 


9480^^ 


T{2S) 


10023 


10004 


10019 


10034 


10024 


10020 


10039^^ 


10025^^ 


r{3S) 


10355 


10350 


10354 


10350 


10349 


10356 


10357 


10352 


Ti4S) 


10580 


10642 


10622 


10585 


10599 


10624 


10593 


10603 


XboilP) 


9860 


9854 


9863 


9854 


9855 


9863 


9868 


9860 


XbiilP) 


9892 


9878 


9898 


9885 


9890 


9900 


9885 


9896 


Xb2{lP) 


9913 


9896 


9921 


9910 


9915 


9924 


9920 


9921 


Xbo{2P) 


10232 


10213 


10228 


10228 


10226 


10230 


10239 


10232 


Xbi{2P) 


10255 


10236 


10243 


10244 


10242 


10245 


10254 


10248 


Xb2{2P) 


10268 


10254 


10255 


10257 


10256 


10258 


10266 


10262 



^States with equal mixtures of scalar and vector couplings. 

^Here's the flavour-dependent case. 

'^Here's the flavour-independent case. 

'^The hyperfine splitting term is not considered here. 



30 



TABLE IV. Salpeter bb quarkonium masses (in MeV) for the Cornell, Song-Lin and Turin 
potentials 



State^ 


[45] 


Cornell 


Error % 


Song-Lin 


Error % 


Turin 


Error % 


l^^o^ 


9398 


9342 


0.60 


9344 


0.57 


9343 


0.58 


13^1 


9460 


9438 


0.23 


9464 


0.04 


9457 


0.03 


l'P2 


9913 


9912 


0.01 


9900 


0.13 


9907 


0.06 


l3Pi 


9891 


9886 


0.05 


9872 


0.19 


9879 


0.12 


I'Po 


9870 


9846 


0.24 


9834 


0.36 


9838 


0.32 


I'Pi 


9908 


9894 


0.14 


9878 


0.30 


9887 


0.21 


l^Ds 


- 


10149 




10162 




10157 




l^F, 


- 


10335 




10356 




10348 




I'G, 


- 


10496 




10514 




10507 




2^ So 


9983 


9986 


0.03 


9983 


0.0 


9984 


0.01 


2^ Si 


10023 


10008 


0.15 


10024 


0.01 


10015 


0.08 


3^ Si 


10355 


10351 


0.04 


10349 


0.06 


10347 


0.08 


4^ Si 


10580 


10623 


0.41 


10588 


0.07 


10601 


0.20 


Fits^ 


A 


0.470 




0.869 GeFi/2 




0.620 GeFi/4 






K 


0.186 GeV^ 




0.558 GeF3/2 




0.304 GeV^/^ 






Vo 


-0.802 GeV 




-0.893 GeV 




-0.823 GeV 






rub 


5.17 GeV 




5.20 GeV 




5.18 GeV 





^States with equal mixtures of scalar and vector couplings. 
'^Recently searched by ALEPH Collaboration [33]. 
'^Here we cite Ref. [10] for the parameter fits. 
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TABLE V. Be meson mass spectrum (in MeV) for the Cornell, Song-Lin and Turin potentials, 
with a flavour-independent case. 

State e = 1 e = 1/2 e = e = 1 e = 1/2 e = e = 1 e = 1/2 e = [45] [50] 

l^S'o 6296.9 6284.9 6272.7 6287.2 6269.3 6250.2 6278.6 6263.3 6247.3 6270 6253 

2^So 6852.5 66846.7 6840.8 6855.0 6848.5 6841.9 6850.0 6843.9 6837.7 6835 6867 

3^So 7220.1 7215.9 7211.7 7190.4 7186.7 7182.8 7200.9 7196.9 7193.0 7193 

4^So 7520.3 7516.9 7513.5 7441.1 7438.5 7435.9 7474.4 7471.4 7468.5 

l^Pi 6724.8 6719.9 6714.9 6726.6 6721.8 6717.0 6725.9 6721.1 6716.2 6734 6717 

2iPi 7100.4 7096.6 7092.7 7091.7 7088.5 7085.3 7095.0 7091.5 7088.1 7126 7113 

1^D2 6979.9 6976.5 6973.0 6992.9 6990.3 6987.7 6988.3 6985.3 6982.3 7077 7001 

2^D2 7298.9 7296.0 7293.0 7277.4 7275.4 7273.4 7287.4 7285.0 7282.5 

l^Pi 6717.5 6706.5 6695.4 6724.9 6710.9 6696.5 6719.7 6707.5 6695.0 6749 6729 

2^Pi 7098.4 7090.1 7081.8 7091.7 7082.7 7073.4 7093.1 7084.6 7076.0 7145 7124 

1^D2 6981.4 6974.0 6966.6 6994.5 6987.2 6979.8 6989.4 6982.2 6974.9 7079 7016 

2^D2 7300.4 7294.0 7287.7 7278.7 7273.1 7267.4 7288.4 7282.6 7276.6 

l^Si 6296.9 6300.8 6304.8 6322.7 6327.6 6332.4 6310.1 6314.4 6318.6 6332 6317 

2^Si 6852.5 6854.4 6856.4 6862.4 6864.5 6866.6 6854.4 6856.3 6858.3 6881 6902 

3^51 7220.1 7221.5 7222.9 7193.5 7194.7 7196.0 7202.5 7203.8 7205.1 7235 

4^Si 7520.3 7521.4 7522.6 7442.8 7443.6 7444.5 7475.2 7476.2 7477.2 

1^P2 6735.9 6748.5 6762.0 6739.2 6755.1 6771.8 6747.0 6753.6 6768.0 6762 6743 

2^P2 7104.7 7115.6 7126.4 7097.9 7109.2 7120.2 7101.4 7112.2 7122.8 7156 7134 

l^Ds 6979.4 6998.1 7016.3 6996.1 7013.6 7030.9 6990.4 7008.3 7025.3 7081 7007 

2^D3 7298.3 7314.4 7331.1 7279.4 7293.3 7307.5 7288.9 7303.6 7318.8 

l^Po 6673.1 6639.2 6602.9 6695.8 6652.0 6602.0 6676.7 6636.7 6592.3 6699 6683 

2^Po 7082.4 7060.7 7038.6 7078.0 7052.4. 7025.2 7075.9 7052.5 7028.2 7091 7088 

l^Di 6977.4 6949.6 6921.2 6988.9 6960.1 6929.5 6983.2 6955.1 6926.0 7072 7008 



2^Di 7297.8 7274.1 7250.2 7274.8 7252.9 7230.0 7284.1 7261.6 7238.5 
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TABLE VI. Flavor-independent cb mass spectrum (in MeV) using the Cornell potential with 
scalar (e = 1), an equal mixture of scalar- vector (e = 1/2) and vector (e = 0) interactions. 



State 


e = P 


e = 1/2^ 


6 = 0'= 


State 


e = 1 


e = 1/2 


e = 


Parameters [25] 


l^^o 


6302.6 


6284.5 


6265.9 


I'Si 


6302.6 


6308.5 


6314.4 


A = 0.437 


2^50 


6884.6 


6875.2 


6865.6 


2^Si 


6884.6 


6887.7 


6890.8 


K = 0.203 GeV^ 


3^50 


7285.1 


7278.2 


7271.2 


33^1 


7285.1 


7287.4 


7289.7 


Vb = Od 


4^50 


7614.6 


7608.9 


7603.1 


4^51 


7614.6 


7616.5 


7618.4 


ruc = 1.321 GeV 


l^Pi 


6740.7 


6732.9 


6725.0 


l^P2 


6746.9 


6773.4 


6796.6 


rub = 4.731 GeV 


2^ Pi 


7152.2 


7146.0 


7139.7 


2^P2 


7155.6 


7175.1 


7196.3 




1^D2 


7019.4 


7013.8 


7008.2 


l^Ds 


7012.2 


7046.2 


7078.8 




2^D2 


7371.0 


7366.1 


7361.2 


2^D3 


7366.5 


7394.9 


7411.0 




l3Pi 


6734.5 


6714.3 


6693.8 


l3Po 


6677.3 


6608.7 


6525.4 




2^Pi 


7151.8 


7136.3 


7120.5 


2'Po 


7132.4 


7091.6 


7049.0 




V'D2 


7024.1 


7010.5 


6996.8 


l^Di 


7024.6 


6973.8 


6921.0 




2^D2 


7375.2 


7363.2 


7351.2 


2^Di 


7376.1 


7332.2 


7286.9 





'^Here e = 1 is the scalar confinement interaction. 

^Here e = 1/2 is the equal mixture of scalar-vector confinement interaction. 
'^Here e = is the vector confinement interaction. 
It's a shifting parameter (in MeV) to be set to fix the the position of M{Bc) and M{B*). 
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TABLE VII. cc and hh mass spectra (in MeV) using the Cornell potential. 



State 


Meson'^ 


Theory [45] 


e = 1 


e= 1/2 


e = 


Meson 


Theory [45] 


e= 1 


e = 1/2 


e = 


l^^o 


r]c 


2979 


3068.1 


3020.8 


2970.7 


Vb 


9400 


9448.3 


9441.3 


9434.2 


2'S^ 


V'c 


3588 


3686.7 


3659.6 


3631.9 


< 


9993 


10021.5 


10018.3 


10015.0 


31^0 


v'^ 


3991 


4135.4 


4115.0 


4094.2 


< 


10328 


10380.4 


10378.1 


10375.9 


41^0 






4510.7 


4493.8 


4476.6 






10669.8 


10668.0 


10666.2 


I'Pi 


he 


3526 


3505.7 


3482.7 


3459.3 


h 


9901 


9901.3 


9898.4 


9895.5 


2^Pi 


K 


3945 


3974.1 


3955.5 


3936.7 


K 


10261 


10266.0 


10263.9 


10261.7 


l^D2 




3811 


3812.7 


3796.0 


3779.1 




10158 


10150.5 


10148.5 


10146.5 


2^D2 






4218.2 


4203.6 


4188.8 






10456.4 


10454.8 


10453.1 


l3Pi 


Xcl 


3510 


3500.6 


3463.7 


3425.7 


Xbi 


9892 


9894.9 


9890.5 


9886.0 


23Pi 


X'cl 


3929 


3975.5 


3946.9 


3917.8 


xL 


10255 


10263.7 


10260.5 


10257.2 


13^2 




3813 


3820.8 


3795.6 


3770.0 




10158 


10150.1 


10147.1 


10144.1 


23^2 






4225.4 


4203.2 


4180.9 






10456.3 


10453.8 


10451.3 


13^1 


J/iJ 


3096 


3068.1 


3083.3 


3098.2 


T 


9460 


9448.3 


9450.6 


9452.9 


23^1 


^' 


3686 


3686.7 


3695.6 


3704.4 


T' 


10023 


10021.5 


10022.6 


10023.7 


33^1 


r 


4088 


4135.4 


4142.1 


4148.8 


-y// 


10355 


10380.4 


10381.2 


10381.9 


43^1 


r 




4510.7 


4516.3 


4521.9 


-y/// 




10669.8 


10670.4 


10671.0 


13P2 


Xc2 


3556 


3517.1 


3566.0 


3613.5 


Xb2 


9913 


9910.3 


9916.7 


9923.0 


23^2 


Xc2 


3972 


3978.8 


4019.3 


4059.0 


Xb2 


10268 


10269.9 


10274.7 


10279.5 


13^3 




3815 


3802.8 


3871.4 


3937.8 




10162 


10153.1 


10161.5 


10169.8 


23^3 






4210.1 


4270.6 


4329.5 






10458.2 


10465.1 


10472.1 


I'Po 


XcO 


3424 


3369.4 


3106.0 


- 


Xbo 


9863 


9864.0 


9850.7 


9837.2 


23Po 


XcO 


3854 


3937.0 


3849.1 


3752.9 


X'bo 


10234 


10251.9 


10243.0 


10234.0 


l^D, 




3798 


3817.8 


3711.0 


3594.8 




10153 


10144.3 


10132.4 


10120.4 


2^Di 






4224.0 


4132.6 


4035.6 






10452.3 


10442.5 


10432.7 



^Same parameter fits as in Table VI. 
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TABLE VIII. Pseudoscalar and vector decay constants (/p = /bc fv = fs*) of the Be meson 
(in MeV). 

Constants SLNET^ rel [45] NR [45] [23] [29] [25] [48] 

fs, 378-513 433 562 479-687 460±60 517 420±13 

/sj 347-482 503 562 479-687 460±60 517 

^We used different potential models with (e = 0). 



TABLE IX. The radial wave function at the origin (in GeV^) calculated in our model and by 

the other authors. 

Level SLNET^ Martin EQ[23] F[25] [23]'" 

\RbM\'^ 0.94-1.73 1.716 1.638 1.81 1.508-3.102 

|2 



Rb*{0)\ 0.73-1.52 



*We used different potential models with (e = 0). 
'"For the IS level. 
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